Magneto-Josephson effects in junctions with Majorana bound states 
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We investigate ID quantum systems that support Majorana bound states at interfaces between 
topologically distinct regions. In particular, we show that there exists a duality between particle- 
hole and spin degrees of freedom in certain spin-orbit-coupled ID platforms such as topological 
insulator edges. This duality results in a spin analogue of previously explored 'fractional Josephson 
effects' — that is, the spin current flowing across a magnetic junction exhibits 4ir periodicity in the 
relative magnetic field angle across the junction. Furthermore, the interplay between the particle- 
hole and spin degrees of freedom results in unconventional magneto- Josephson effects, such that the 
Josephson current is a function of the magnetic field orientation with periodicity 4n. 



The possibility of observing Majorana zero-modes in 
condensed matter has captured a great deal of attention 
in recent years. Much effort in this pursuit presently fo- 
cuses on spin-orbit-coupled ID wires, which are closely 
related to edges of 2D topological insulators (TIs). In 
either setting Majorana modes are predicted to localize 
through the competition between superconducting prox- 
imity effects and Zeeman splitting [THE]. Remarkably, 
zero-bias conductance anomalies [THTT] possibly origi- 
nating from Majorana modes have even been measured 
[121 113] very recently in quantum wires. Numerous other 
fascinating phenomena tied to Majorana fermions have 
also been explored, including non-Abelian statistics [TH - 
fTB] . electron teleporation [IT], and exotic Josephson ef- 
fects [Dm EH]. 

Particularly interesting to us here are the Majorana- 
related Josephson effects in quantum wires and TI 
edges. Consider two Majorana modes hybridized across a 
Josephson junction formed by topological superconduct- 
ing regions separated by a narrow barrier as shown in 
Fig. [ljb). The energy splitting of these Majoranas de- 
pends periodically on half the phase difference between 
the right and left superconductors, (0 r — <fii)/2, giving 
rise to a Josephson current with 47r periodicity in <j) r — <f>i 
[D EE] . If, in addition, a third superconductor contacts 
the middle domain, a difference between its phase and the 
average phase {<p r + 4'i)/2 induces a non-local three-leg 
"zipper" Josephson current that divides equally between 
the two leads and is also Aw periodic in <f> r and 4>i [?]■ 
These 'fractional Josephson effects' provide smoking-gun 
signatures of Majorana modes. 

Our claim is that physical quantities of Majorana junc- 
tions in wires and TI edges can also possess 4-7r-periodic 
dependence on the orientations of Zeeman fields applied 
in the plane normal to the spin orbit direction. Notably, 
in some domain configurations the Majorana-mediated 
Josephson current reverses sign after a full 2-7T rotation 



of the magnetic field orientation on one side of the junc- 
tion. Only an additional 2-7T rotation restores the currents 
to their original direction. Thus the mixing between the 
particle-hole and spin degrees of freedom leads to an un- 
conventional magneto- Josephson effect through the cou- 
pling of Majoranas. 

Additionally, the Majorana modes produce a 'spin 
Josephson current' between the magnets providing the 
Zeeman energy, which could also be An periodic in the 
field orientations. Define 9 S as the angle between the wire 
and the Zeeman field at domain s. Spin Josephson cur- 
rents, j , are equivalent to torques (driven partly by the 
Majoranas) that the wire domains apply on the external 
magnets [2 8) . Therefore, they are given by the derivative 
of the system's energy with respect to the magnetic field 
orientations 0: 

f = (1) 

with H being the system's Hamiltonian. In the case of 
TI edges, the spin currents arise as the exact duals of 
Josephson currents, and the orientation of the B-field is 
the exact dual to the superconducting phase (indeed, the 
Josephson current is given by = jr ) [2"9] . We 
emphasize that the 4tt periodicity prevails as long as the 
parity of the Majorana state remains constant during the 
measurement, or changes at a slower rate than the wind- 
ing of the superconducting phase and magnetic orienta- 
tions. 

Let us focus first on the analysis of the 4-7r-periodic ori- 
entation dependence in TI edges, before commenting on 
spin-orbit-coupled wires which obey qualitatively simi- 
lar rules. The Hamiltonian, including s-wave pairing and 
Zeeman fields in both the transverse and parallel direc- 
tions relative to the spin-orbit direction, reads 

% = vpr z a z - [it z + A (cos <f> t x - sin <f> t v ) 

-ba z +B (cos 9 a x -sin (9 a v ) . (2) 
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FIG. 1: (a) Phase diagram for ID system: gapless- 
phase (B 2 < and A 2 < b 2 ), A-phase (A 2 - b 2 > 
max [B 2 - n 2 , 0] ), and B-phase (B 2 -/i 2 > max [A 2 - b 2 , 0] ). 
Both A-phase and B-phase are gapped, (b) The A-B-A junc- 
tion supports Majorana bound states at the domain walls [4] . 
(c) The dual configuration of B-A-B junction that also sup- 
ports Majoranas. 



Here we have employed the Nambu spinor basis ~^> T — 
(tpf, ipi, tpp — i/>j) and introduced Pauli matrices a a and 
T a that act in the spin and particle-hole sectors, respec- 
tively. The edge-state velocity is given by v, p is the mo- 
mentum, and the (redirection represents the spin-orbit- 
coupling axis. We allow the chemical potential /i, su- 
perconducting pairing Ae 1 ^, longitudinal magnetic field 
strength b, transverse magnetic field strength B, and the 
transverse-field orientation angle 9 to vary spatially In- 
terestingly, Eq. ([2]) has a magnetism-superconductivity 
duality — the Hamiltonian takes the same form upon in- 
terchanging the magnetic terms {b, B. 8, a a } with the su- 
perconducting terms {/u, A, tfi, r a }. Below we deduce the 
physical consequences of this duality. 

The Hamiltonian ^ supports three different phases 
determined by the relative strength of {A, /z, B,b}. As 
Fig. [T|a) illustrates, we have (i) a topological supercon- 
ducting gapped phase (denoted henceforth as the A- 
phase) when A 2 — b 2 > max [B 2 — p? , 0] , (ii) a topo- 
logical magnetic gapped phase (denoted B-phase) when 
B 2 — (i 2 > max [A 2 — 6 2 ,0], and (iii) a trivial gapless 
state when B 2 < /.i 2 and A 2 < b 2 . Consistent with the 
magnetic-superconducting duality, in the phase diagram 
of Fig. [ija) the B- and A-phases are symmetrically ar- 
ranged with respect to the diagonal line that defines the 
boundary between these two gapped states: 



B 2 + b 2 . 



(3) 



Majorana zero-modes bind to domain walls separating 
B- and A-domains. For notational simplicity, below we 
will assume that A > b > and B > \l > 0, though 
more general results can be obtained [15]. We will also 
focus on setups for which all domains experience both 
superconductivity and a transverse Zeeman field. 

In TI edges, the 4tt periodic dependence on the mag- 
netic field orientation occurs when two Majoranas are 
nestled in a B — A — B domain sequence as in Fig. [ljc). 
This is in contrast to the previously studied uncon- 
ventional Josephson effects [TJ UJ [TB] , which occur over 



a junction between two A-domains bridged by a B- 
domain [see Fig.[IJb)]. The magneto- Josephson and spin- 
Josephson effects of a TI edge follow from the detailed 
dependence of the Majorana energy splitting, i?Maj ; on 
the field orientations and superconducting phases in the 
B — A — B edge domain structure of Fig. [IJc) . In addi- 
tion to an exact numerical calculation of i^Maj , we pro- 
vide in |19j an analytical variational approach that sheds 
light on the physics. In the latter approach we assume 
that the Majorana wavefunctions are unmodified by their 
proximity to each other, apart from being superposed to 
form a conventional low-lying state. This leads to an en- 
ergy splitting that is suppressed as a weighted sum of 
two exponentials which control the decay of the Majo- 
rana wave functions in the middle domain. 

Our results for the Majorana couplings constitute one 
of the central results of this paper. The two characteristic 
decay lengths as a function of field and pairing are Ai^ = 
— -, Quite generally, for the middle A- 
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domain of length L, the Majorana coupling energy is: 
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s> 2- K /A 2 rt -b 2 n (A 2 n +lJ, 2 n -B 2 n -b 2 rl ) 
(Bs+bs-Ms) + A s [y/B 2 -^ 2 s\n{l m ±S(j> s )-b s cos{b m ±Scj> s )] 

with the choice of sign ± depending on s = I or r. Note 
that M s exhibits the standard 2-7T periodicity in (f> s , so 
that the more exotic 4n periodicity follows exclusively 
from the trigonometric functions in Eq. Q. 

These general results allow us to quantitatively es- 
timate the magneto- Josephson effects described ear- 
lier, which can be measured in the circuit sketched in 
Fig. [2]J a). For simplicity, we specialize to the case of 
f-i/m/r — 0, where the Majorana coupling energy reduces 
to 
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The Major ana-related magneto- Josephson currents en- 
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where p = ±1 denotes the parity of the hybridized Ma- 
joranas. The explicit form for the charge currents (drop- 
ping the parity factor p) is: 
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which constitutes a prediction for the unconventional 
magneto- Josephson effect. The analytical expressions ob- 
tained above for jfj and j% agree well with the numeri- 
cal calculations for large L as shown in Fig. [2jb) . They 
confirm that for B-A-B junctions the Majorana coupling 
induces the charge current jy r with An periodic depen- 
dence on 9i/ r . 

Similarly the spin Josephson currents, or torques on 
tne magnets, m region s — l/r are j£ = — = p dg ' 
[Eq. Q]. The angular momentum transferred by these 
currents is in the direction parallel to the spin-orbit axis, 
which in this case is the z-direction. The spin Josephson 
currents are thus given by: 



aS i_„-S „;„ 9i — 8 r , -S „:„ 8i+9,, — 29„ 

3l/r ~ ^3m Sm 2 +3z sin 2 

with: f M/z = ^. 



(9) 



The jff spin current exchanges angular momentum be- 
tween the right and left magnets directly, while the j§ 
spin current originates in the middle region and equally 
splits into the right and left regions, j^ n ^ t — 3m^ r = 
jS gm 9i+9 r -29 m ^ This t erm vanishes when there is no 
transverse magnetic field in the middle domain, and rep- 
resents the dual of the zipper Josephson effect in the A- 
B-A junction that splits charge current from the middle 
domain between the two side domains [J]. 

The origin of this exotic dependence of the Majorana- 
related currents can be traced to the magnetic- 
superconducting duality in topological insulator edges 
PJ 0]. For a junction with three alternating domains, 
there are two dual configurations: the A-B-A junction 
[Fig. [ijb)] and the B-A-B junction [Fig. [TJc)]. The 
spin- Josephson effect in the B-A-B junction is dual to 
the charge- Josephson effect in the A-B-A junction [T]- 
|4] . Similarly, the magneto- Josephson effect depending on 
the orientation angles in the B-A-B junction has a dual 
spin- Josephson effect depending on the superconducting 
angles in the A-B-A junction. 

Majorana junctions in spin-orbit coupled wires exhibit 
the same magneto- Josephson and spin- Josephson effects 
as the TI edge. The wire's Hamiltonian adds a kinetic 
energy piece to Eq. (|2j), Hk = ^V 2 ' 1 ' 7 ' ' ■ This produces 
additional Fermi points at 'large' momenta pp ~ ±2mv 
that are, however, nearly unaffected by the magnetic field 
in the presence of pairing. Therefore the analysis above 
for the TI edges still applies qualitatively. Thus, in a 
Majorana wire, 47r-periodic effects in both 9 and 4> appear 
in the B — A — B domain sequence [3D] . The quantitative 
analysis of the magneto-, spin-, and charge- Josephson 



TI edge 



SC 







, — — > 




, — Q > 


> 


























. > 

)m 




t . 

Jz 


1 


H 




Jz 




j S r (t) = yb r t 

Jm 




FIG. 2: (a) The scheme to measure unconventional magneto- 
Josephson effect. Josephson currents are measured for the 
B-A-B junction. In the right region, the transverse magnetic 
field winds at rate ljl = yb r , which modulates the Josephson 
current at half the frequency, Wl/2. (b) Comparison between 
analytical expressions and numerical results for jm and jz- 
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effects in wires as well as the role of Andreev bound states 
will be analyzed elsewhere pp] . 

Observing the unconventional magneto- Josephson ef- 
fect and the Air periodicity in 8 t / r [see Fig. p^b)] requires 
effective control of the magnetic field orientation. In par- 
ticular, the orientation change needs to be sufficiently 
fast so that the Majorana states' total parity does not 
change, but still slow on the scale of the inverse bulk 
gap to avoid quasiparticle poisoning [3TJ. The rate of 
parity decay is strongly detail dependent, but we sur- 
mise that measurements with rates faster than 1 kHz and 
slower than the minimum gap in the device would suffice. 
Conventional magnets may be too unwieldy when made 
to rapidly turn; nuclear magnetization, however, could 
be ideal for this task. Through the hyperfine coupling, 
a polarized nuclear spin population could create an ef- 
fective Zeeman field in the plane perpendicular to the 
spin-orbit coupling direction. For example, large nuclear 
spin polarization, normal to the spin-orbit direction, can 
be induced by optical pumping with circularly polarized 
light [35]. An external magnetic field with strength 6, 
applied parallel to the spin-orbit axis, would make the 
orientation angle of the hyperfine transverse field wind 
at a rate uj l = 76, where j/2tt « -7.6 MHz/T for 199 Hg 
or 7/2tt rj 13.5 MHz/T for 125 Te nuclei [53]. The hyper- 
fine transverse field can be rather strong, e.g., B ~ 0.1 
Tesla for 2% nuclear polarization fraction [22]. It can, 
moreover, persist for long times, limited by the inhomo- 
geneous nuclear transverse spin lifetime ~ 100 /is, 
which already suffices for hundreds of precession periods 
for b ~ 0.1 Tesla. The transverse spin lifetime can be 
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FIG. 3: Contour plot of (a) spin current jr and (b) charge 
current , both of which are 4-7T periodic in 9 r and 2-zr periodic 
in cj> r . The other angles are fixed <f>i = n/2, cj> m = 6i/ m = 0. 
The parameters are the same as Fig. [2] 

further extended using spin echo techniques. 

With a rotating transverse magnetic field, we can 
observe the magneto- Josephson effect in several ways. 
A constantly winding orientation in the left domain, 
6 r (t) = oj^t [while fixing 9[/ m (t) = 0, as illustrated 
in Fig. [2^a)], produces an oscillatory component of the 
charge current with amplitude j® L / 2 — 3 m + 3z ~ 

TlW e ~ LIXm ~ at hal f thc frequency, u L /2. In TI 
edges, we can also use resonant properties to probe the 
orientation-frequency halving. A DC voltage V applied 
to the right superconducting lead, for instance, induces a 
winding of the superconducting angles, <p r (t) = 2eVt/h 
and (f>i/ m (t) = 0. When the magnetic orientation also 
winds with angular velocity f2^, interference between the 
two oscillation would yield a DC current from the right 
superconducting lead, when cjj, = 2Qy (neglecting high- 
order resonances). The amplitude of the dc current is 
expected to be: 

1 /" 2lr 

3Zf=2n v ~ ^ y o J^/ 2 [<Mcos<MaV. (10) 

Alternatively, one can apply an AC voltage to the right 
superconducting lead such that <f> r oc sinwi, while all 
other superconducting angles are held fixed. Interference 
effects now produce Shapiro-step-like resonant features 
which emerge only when 

Wi = 2nw (11) 

for even integer 2n (neglecting higher order corrections 
to the dependence). 

The Majorana-mediated spin currents with 47r phase 
periodicity are harder to measure. A possible route for 



such measurements is to use a magnetic nanoparticle as 
the magnetic field source on one of the side domains. The 
torques on the nanoparticle could be probed from the 
shift in the ferromagnetic resonance (FMR) frequency. 
The FMR frequency is typically f ~ 10 GHz. The FMR 
linewidth, dictated by the Gilbert damping coefficient a, 
is of order cv/o = 0.01/o in bulk ferromagnets, but is 
probably much smaller in nanoparticles 24 . A rough 
estimate of the maximum Majorana-related spin-current 
(or torque), j s , yields j s ~ H- 10 GHz. This produces a 
frequency shift around j /mtotai, which is inversely pro- 
portional to the total angular momentum of the FM grain 
mtotai [25]. This shift must dominate the FMR linewidth, 
j /mtotai > foot- The nanogram must, therefore, be suf- 
ficiently small such that m tota i/H < a^ 1 ~ 100, e.g. have 
a radius of around 10 nm, and still provide a sufficient 
Zeeman field for the domain it is on. 

Measuring the effect of the relative field orientation on 
the spin and charge currents can be complicated by the 
presence of conventional Josephson effects arising from 
the continuum states. Indeed, the bulk energy associ- 
ated with the continuum states also has dependence on 
magnetic field orientations and superconducting phases 
that are interesting in their own right, and of similar mag- 
nitude to the Majorana related effects. Nonetheless, all 
these dependencies are 2-7r periodic, as we have confirmed 
numerically. Hence, the measurement schemes proposed 
above will be insensitive to them. 

In conclusion, we explored consequences of a 
magnetism-superconductivity duality of TI edge states, 
emphasizing Josephson effects. Most prominently, the 
duality implies that spin and charge Josephson currents 
in TI edges exhibit a 4-7T periodic dependence on the orien- 
tation difference of the magnetic field. These remarkable 
effects are a direct consequence of the Majorana states 
and we make several proposals how to detect them exper- 
imentally. The duality is only approximate in spin-orbit- 
couplcd quantum wires but analogous effects also occur in 
this system. In addition to the Josephson effects, the du- 
ality has further interesting implications. For instance, it 
implies that the transition between topological and triv- 
ial phases can be tuned using a magnetic gradient, which 
is the dual of the superconducting phase gradient 26 . 

Note added: As we are completing the manuscript 
we became aware of overlap work by Qinglei Meng, Va- 
sudha Shivamoggi, Taylor Hughes, Matthew Gilbert, and 
Smitha Vishveshwara [27] . 
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perconductor. If the p = crossing forms another p-wave 
superconductor, together the two form a topologically 
trivial phase. 
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Supplementary materials 

We study the linearized ID system 

H (fi, A, <j);b,B,e)= pr z a z - fir z + A (t x cos - t v sin 0) - ba z + B (a x cos 9 - <r y sin 0) (12) 

with six control parameters: fi for the chemical potential, A for the pairing energy, <f> for the superconducting phase, 
—b for the longitudinal magnetic field, B for the transversal magnetic field, and 6 for angle of the transversal magnetic 
field. In this form, the duality between (A, ft) and (B, b) is more obvious. Without loss of generality, we assume that 
all the control parameters (//, A, b, B) are all positive. 



Phase Diagram 



We compute the determinant 



detH 



P 



Vb^ 



Va 2 - b 2 



p 



Vb~i~ 



a/A 2 - b 2 



(13) 



The energy gap will be closed if there exist some real solutions of p to satisfy det H = 0. 

1. When B 2 < fi 2 and A 2 < b 2 , the system is in a gapless-phase, because there are real solutions p — 
± [y/-B 2 + p? + \/-A 2 + 6 2 ) or p = ± [y/-B 2 + p 2 - \/-A 2 + b 2 ^j to fulfill the requirement of det H = 0. 

2. When B 2 > p 2 or A 2 > b 2 , the system is always gapped, because there are no real solutions of p to satisfy 
det H = 0. 

(a) For A 2 — b 2 > max [B 2 — /i 2 , 0] , the system is in a superconducting gapped phase (A-phase). 

(b) For B 2 — p? > max [A 2 — b 2 , 0] , the system is in a magnetic gapped phase (_B-phase) . 

(c) There is a quantum phase transition at A 2 — b 2 = B 2 — fi 2 , which connects the A-phase and the J5-phase. 

Therefore, we obtain the phase diagram in Fig. [lja). 



ID System Consisting of Different Regions 



We are interested in the case that the ID system consists of three regions of different control parameters. Specifically 

{Xi for x £ (— oo, 0) 
Xm forze (0,L) ' (14) 
Xr for x € (L, +oo) 

with x representing the six control parameters. The system Hamiltonian is 

(Hi for x e (-oo, 0) 
H m for x £ (0, L) (15) 
H r for a; G (L, +oo) 

with Ht = H (pf,Af,(f>f-,bf,Bf,6f). We are interested in the B — A — B configuration, with Bf — fif > 
max [A 2 - bf, 0], A 2 n -b 2 n > max [B^ - p? m , 0] , and B 2 r - p? r > max [A 2 - 6 2 , 0] . 



Perturbative Formulism for the Coupling Energy 



Let's first consider the individual Majoranas. The left Majorana \L) is at x — associated with the l — m boundary. 
We may introduce the Hamiltonian Hl = I ^ X ^ / n 00 l^' that supports the zero energy Majorana mode \L), 
with Hl \L) = 0. Similarly, the right Majorana \R) is at x = L associated with the m — r boundary. We can also 
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{H for x £ ( oo T j i 
™ f )t i \ that su PP or ts zero-energy Majorana mode \R), with Hr\R) = 0. We can 
rL r tor x £ (-^j) too J 

can perturbatively compute the coupling energy between \L) and \R) by the formula: 

Hlr « M" 1/2 /iAr 1/2 (16) 
with M being the overlap matrix between the (not necessarily normalized) Majorana states, and h being: 



h = 



(L\AV\R) 
(R\AV\L) 



with: 



Therefore, the coupling Hamiltonian is approximately Hlr ~ ( ^* q ] , with 



AV = H - H L = (H r - H m ) n{x-L). 

£ 



(Lj AV L \R) 
y/{L\L) (R\R) 



(17) 



(18) 



(19) 



Wavefunction of Individual Majoranas 



We can rewrite the Hamiltonian as 



H L = H lV (-x) + H mV (x) 

U r V-(p-Ki)-V^ -U\t z u z fora:<0 
U m ■ V ■ (p - K m ) • W • Ulr z a z for x > 



where the unitary transformations are 



V = e 



(20) 
(21) 

(22) 



U = e^ T ® = : Uj, <g> C/ e 

and the non-Hermitian matrix is 

A" = (6t z + iAt x ) + (fia z + iBo x ) 
with sub-index / = l,m,r not explicitly written for simplicity. Without loss of generality, we can fix 





(23) 
(24) 
(25) 



and U m = I. For our notational convenience, we also introduce b = cos 1 ^ and ft = cos 1 -g. (Let's assume A 2 > b 2 
and B 2 > /i 2 for notational simplicity. Later we will show that this constraint can be relaxed.) The eigensystem of 
K is 



K ■ 



with sub-eigenvectors 



(26) 



(27) 



and eigenvalues 



A SUS2 = A A? 1 + B A~ 2 



(28) 
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where 



A| = = i sin s£ 



for s = ±1. ■ v~£; = —i sin ^iS-. The two-vectors w| = have the following properties of inner-products: 



(w|) T • u| = — isins£ <5 S; ,. 



(t>|) t • a z ■ = -isins£ 5s, s < 



vf = 5 ata i + COs£ 5s, s 



-isin£ 

i sin £ 

— isin£ 
i sin £ 



1 cos£ 
cos£ 1 



where s := — s for s = ±1. And it transforms under the unitary 



Uevj 



Left Majorana. 

For the B — A interface at x = 0, the localized zero-energy eigenstate is 

Ui ■ r z a z |* Q ) fora;<0 



\V -U m -T z a z \^p) forx>0 ' 



with 



One can verify 



because 



s 
s 

#L |i> = 



f (p-if,)l*a>=( 



= for x < 
for x > 



The boundary condition \L(x = )) = \L (x = + )) requires 

L/i|* a (S = 0-)) - (z = + )), 

and hence 

J2 v l a s = u-^ 



which gives us 



a s = sin 1 sfe; sin 



a . -1 - . SMm + (Ol+fk) 
p s = sin s^ m sin . 



(29) 

(30) 
(31) 
(32) 

(33) 



(34) 

(35) 
(36) 

(37) 
(38) 

(39) 
(40) 
(41) 

(42) 
(43) 
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Right Majorana. 

Similarly, For the A — B interface at x = L, the localized zero-energy eigenstate is 



with 



and 



, . _ f V ■ U m ■ r z a z |* 7 ) for x < L 
1 ; \ V ■ U r ■ t z g z |* 4 ) forx>i? 

* 7 (x) = ® v' a 7 a e- iA &'(*- L ) 

s 



. _1 - . sAm + {Or - Ar) 

7 S = sin S/U m sin 



5 S = sin 1 sbi sin ■ 



s6 r - (j) r - b rn ^j 



(44) 

(45) 
(46) 

(47) 
(48) 



Normalization of Wavefunctions 



The normalization of wavefunction is 
M l = (L\L) 



poo 

= dx(*f}(x)\*(}(x))+ / dx{9 a (x)\9 a {x)) 

JO J -co 

~ 2 VA? n - fe2„ (A*, + ^ - ift - ft*,) . 

{Bf + bf- tf) + A t [y/Bf - sin (b m + </>,) - h cos (b m + 0,) 



+ 



2^Bj^{Bf + bf-Af-tf) 



(49) 



Note that the each of the two terms are positive definite, because Bf + bf > Af + [if and A^ + ^„ > B^ + b 2 m . By 
taking [if — (i.e., /ty = ir/2), bi^ m ^ r = (i.e., 6; im>r = 7r/2), we have the expressions 



<i|i>c 



+ 



Br + A; COS ( 



/0 ° 2(A^-i^) ' 2(B?-A*) ' 
We can also compute M r [<j> r ] = {R\R}, which is very similar to Mi [</>;] with the following replacement 

b m + (j>i =>■ b m - <j) r 
Ai,m,Bi,bi =>■ A r , fi r ,B r , b r 



(50) 



(51) 
(52) 



Cross Coupling (L\ AV L \R) 

We now compute the cross coupling term (L\ AVl \R). First, we can rewrite AVl 

AV L = -U r ■ V ■ {p - K r ) ■ V* ■ U\t z g z xrj{x-L) 
+ U m -V-{p-K m )-V^Ul- t z o z x V {x-L). 



(53) 



10 



The matrix element 
(L\ AV L \R) 

dx (9 p (x)\ Ul-U r - T z a z ■ K r |¥ 4 (x)) + \ dx (K m )* ■ r*a z ■ W m ■ U r |tf 4 (x)> 

2 <)E^- iA ^*^> 2 



sm /i,, 



zMmzzMz 



(54) 



By taking r = (i.e., /2;, m ,r = 7r /2) ) fy,m,r = (i.e., &j, m , r = tt/2), we restore the previously obtained familiar 
expressions 



(L| AVi |i?) 00 ex e 



-A m L / BmL "I "r , — B™ L ■ °l ■ u r 
2 2 2 2 



(55) 



Majorana Coupling Energy 

We can compare the perturbative calculation with the numerical results. For simplicity, we choose the parameters 
fJ-l/m/r = 0) h/m/r = 1/2, A m = 2.5, A;/ r = 1, -E>;/ r = 2, £? m = 1. The energy from perturbative calculation is 



E 



Maj 



(L\ AV L \R) 



a/M; [^] M r [<jV] sin //„ 



oVSV^i sin 9 -+A 2 "+* sin g^fe^ 



+e -V B ™-^ L sin ft-A^+A. sil 



2 

Am -Mr- 



(56) 



For this set of parameters, it will be better to choose 0/ = 7r/2, 4> r — tt, so that E will be most sensitive to the 
deviation in <p, which gives the max charge current Iq cx ^M. 



Analytic Continuation for A 2 < b 2 or B 2 < fi 2 



For A 2 < b 2 or B 2 < fi 2 , we may define the complex number from the analytic continuation 

I -1 b u-i 5 

6 = cos — = — z cosh — 

A A 



(57) 



-i /•* - , -i I* 
ix = cos — = — 2 cosh — . 

B B 



(58) 



The eigensystem of K has the same form K ■ Ov- 1 ® w~ 2 ^ = A Sl,S2 fv? 1 ® «? 2 \ with sub-eigenvectors 



^2 



£/ 2 , e IS £/ 2 ) = u^, eigenvalues A S1,S2 = A A- 1 + B A~ 2 , and sub-eigenvalues 



A| = A^ = i sin s£ = sinh sz£ 

for s = ±1. The orthogonality condition remains consistent with the analytic continuation: 

/ S \T s > . . . r / — sinh if 

(« e ) » Bin ( o sinhif 



(59) 



(60) 



Hence, the coefficients {a s , /3 S , -f s , 5 S } can be obtained by analytic continuation from Eqs.( 42|43|47|48 ). For example, 



a s = sinh 1 sibi sinh — - — - . In order to obtain the the overlap of wavefunctions, we need to compute (vf\ . 
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When cos£ > 1 (< 1), e %s( -l 2 is real (imaginary) and (v'^j — — v^-a z (= i v^-a x ). Hence, for cos£ > 1, ■cr z -u| = 

isins£ <5 S)S / and (y^j ■ v| = coshj£ 5 S;S / + <5s,s'- After some careful calculation, we can verify that the analytic 

continuation from Eqs.( 49|54 1 also give the correct results for A 2 < b 2 and/or B 2 < fi 2 . Therefore, Eq.(56l and its 
analytic continuations give the coupling energy between two Majorana bound states. 



